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<N . ABSTRACT 

Many different families of mass models are used in modern applications of strong 
gravitational lensing. I review a wide range of popular models, with two points of 
£NJ ■ emphasis: (1) a discussion of strategies for building models suited to a particular lensing 

problem; and (2) a summary of technical results for a canonical set of models. All of the 
models reviewed here are included in publicly- available lensing software called gravlens. 
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1. Introduction 

As applications of strong gravitational lensing have become more sophisticated, the variety 
and complexity of mass distributions used for lensing studies have increased. Gone are the days 
when the singular isothermal sphere was all you needed to know. Now we have softened power law 
ellipsoids, pseudo-Jaffe models, generalized NFW models, exponential disks embedded in isothermal 
halos, and so on. Plus, all of these models are likely to be perturbed by other galaxies, groups, or 
clusters near the lens galaxy or along the line of sight. 
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I have developed software called gravlens for a variety of lensing and lens modeling applica- 
tions. 2 In the course of writing the code I have collected most of the mass models used for lensing 
studies, and this paper reviews those models. The outline is as follows. Section 2 discusses circular 
and elliptical symmetry in lensing mass distributions, and argues that a canonical set of circular 
and elliptical models provides a useful basis set for building much more complex composite mod- 
els. Section 3 gives some suggestions (propaganda, really) about selecting models appropriate for 
your application. Section 4 presents the general equations that describe the lensing properties of 
a given mass distribution. Finally, Section 5 presents a detailed catalog of results for a number of 
models. The mass models discussed in this paper (and included in the gravlens software) are listed 
in Table 1. 
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2 The software is discussed in a separate paper (Keeton 2001), and is available to the community via the web site 
of the CfA/Arizona Space Telescope Lens Survey, at http://cfa-www.harvard.edu/castles. 
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2. Circular, Elliptical, and Composite Lens Models 

The lensing properties of any mass distribution can be written in terms of two-dimensional 
integrals over the surface mass density (see §4). In general the integrals cannot be evaluated 
analytically, but many lensing applications offer simplifications due to symmetry. Applications 
such as microlensing may permit the use of a lens with circular symmetry, in which case the lensing 
properties can usually be found analytically In other applications, such as lensing by galaxies, it 
may be reasonable to assume elliptical symmetry, which allows the lensing properties to be written 
as a set of one-dimensional integrals (see §4); the integrals can sometimes be evaluated analytically 
and are always amenable to fast numerical techniques. 

The geometric simplifications may not hold in increasingly sophisticated lensing applications. 
For example, companion stars or planets can break circular symmetry in microlensing (e.g., Mao 
& Paczyhski 1991; Gould & Loeb 1992); neighboring galaxies (e.g., Young et al. 1981; Schechter et 
al. 1997; Koopmans & Fassnacht 1999) or the lens galaxy's internal structure (e.g., Mailer, Flores 
&; Primack 1997; Bernstein & Fischer 1999; Keeton et al. 2000) may break elliptical symmetry in 
lensing by galaxies; and individual cluster members can ruin the symmetry of the smooth back- 
ground potential in lensing by clusters (e.g., Bartelmann & Steinmetz 1995; Tyson, Kochanski h 
Dell' Antonio 1998). Nevertheless, in all of these examples the total mass distribution can be written 
as a combination of circular or elliptical components, either placed at different positions to repre- 
sent different objects, or combined at the same position to mimic an object with more complicated 
internal structure (also see Schneider & Weiss 1991). Because the Poisson equation is linear, the 
lensing potential of such composite models is simply the sum of the component potentials. In other 
words, arbitrary combinations of circular and elliptical models yield a wide range of complex mass 
distributions whose lensing properties are nevertheless easy to compute. 

While composite models provide a great deal of freedom and complexity, they are not com- 
pletely general. This limitation is eliminated in the elegant algorithm by Saha & Williams (1997; 
also Williams & Saha 2000) for finding non-parametric lens models. The approach is to introduce a 
set of mass pixels and construct a large linear programming problem for determining their masses. 
The problem is severely underconstrained, but by requiring positive-definite masses and imposing 
some smoothness criteria it is possible to find a wide but finite range of non-parametric models 
consistent with the data. 

A limitation of the Saha & Williams algorithm is that the constraint of a positive-definite 
surface density is weaker than the constraint of a positive 3-d mass density — or, better yet, a 
positive-definite quasi-equilibrium distribution function. In other words, while parametric models 
may provide too little freedom, the Saha & Williams non-parametric models provide too much. For 
example, their method does not determine whether models are consistent with stellar dynamics, 
and many of the models found by the method probably are not consistent. The gravlens software 
includes only parametric models that have physical motivations apart from lensing, but it offers 
the ability to combine them in arbitrary ways to achieve extensive and reasonable complexity. 
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3. Selecting Models 

Selecting the class of models to use for a particular application is a key part of lens modeling. 
When modeling lenses produced by galaxies, a simple and useful place to start is an isothermal 
model, i.e. a model with density p oc r~ 2 and a flat rotation curve. Spiral galaxy rotation curves 
(e.g., Rubin, Ford & Thonnard 1978, 1980), stellar dynamics of elliptical galaxies (e.g., Rix et al. 
1997), X-ray halos of elliptical galaxies (Fabbiano 1989), models of some individual lenses (e.g., 
Kochanek 1995; Cohn et al. 2001), joint lensing and dynamical analyses (e.g., Treu & Koopmans 
2002; Koopmans & Treu 2002), and lens statistics (e.g., Maoz & Rix 1993; Kochanek 1993, 1996) 
are all consistent with roughly isothermal profiles. 

However, an isolated isothermal ellipsoid rarely yields a good quantitative fit to observed 
lenses (e.g., Keeton et al. 1997, 1998; Witt & Mao 1997). In general, adding parameters to the 
radial profile of the galaxy fails to produce a good fit, but adding parameters to the angular 
structure of the potential dramatically improves the fit (e.g., Keeton & Kochanek 1997; Keeton et 
al. 1997). The additional angular structure comes from the tidal perturbations of objects near the 
main lens galaxy or along the line of sight. In other words, the fact that few galaxies are truly 
isolated means that lens models generically require two independent sources of angular structure: 
an ellipsoidal galaxy plus external perturbations. The combination of angular terms can make 
it difficult to disentangle the shape of the galaxy and the nature of the external perturbations, 
and it is extremely important to understand any degeneracies between the two sources of angular 
structure before drawing conclusions from the models (see Keeton et al. 1997). 

To move beyond isothermal models and explore other radial profiles, softened power law lens 
models have traditionally been very popular. However, these models have flat cores, while many 
early- type galaxies have cuspy luminosity distributions (e.g., Faber et al. 1997), and dark matter 
halos in cosmological simulations have cuspy mass distributions (e.g., Navarro, Frenk & White 
1996, 1997; Moore et al. 1998, 1999). The lack of central or "odd" images in most observed galaxy 
lenses also limits the extent to which galaxies can have flat cores (e.g., Wallington & Narayan 1993; 
Rusin & Ma 2001). Cohn et al. (2001) thus argue that softened power law models are outdated and 
should be replaced with families of cuspy lens models. There are three traditional families of models 
in which the profile of the cusp is fixed: NFW (Navarro, Frenk &: White 1996, 1997), Hernquist 
(1990), and de Vaucouleurs (1948) models. In addition, there are three families of models in which 
the cusp is taken to be an arbitrary power law: Nuker law models (Lauer et al. 1995; Byun 1996; 
Keeton 2002), and two variants of NFW models (see §5 for definitions; Jing & Suto 2000; Keeton 
& Madau 2001; Muhoz, Kochanek & Keeton 2001; Wyithe, Turner k Spergel 2001). 

An important reason to study families of cuspy lens models is to test the prediction from 
the Cold Dark Matter paradigm that halos at a wide range of masses are consistent with a unified 
family of halo models (e.g., Navarro et al. 1996, 1997). The distribution of image separations among 
the known lenses rules out the hypothesis that galaxies and clusters can be described by identical 
lens models (Keeton 1998; Porciani & Madau 2000; Kochanek & White 2001). However, it may be 
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possible to resolve the conflict with a model in which galaxies and clusters start out with similar 
halo profiles, but baryonic processes such as cooling modify the inner profiles of galaxies. Kochanek 
& White (2001) show that such a model can match the observed image separation distribution, 
but it remains to be seen whether the model agrees with the lack of central or "odd" images in 
observed lenses, or with detailed models of individual lenses. 



4. General Equations 



For a mass distribution with surface mass density rc(x) = S(x)/S cr in units of the critical 
surface density for lensing, the two-dimensional lensing potential is (e.g., Schneider, Ehlers & Falco 
1992) 

<^( x ) = ^J ln \ x ~y\ K (y) d y ■ W 

The other lensing properties can be derived from the potential. The deflection angle V<p determines 
the positions of images via the lens equation, 



u = x — V0(x) , 

where u is the source position. The magnification tensor, 



du 



l 



4>xx $xy 
i>xy 1 — §yy 



(2) 



(3) 



determines the distortions and brightnesses of images. (Subscripts denote partial differentiation, 
4>ij = d 2 4>/dxidxj.) Many lensing applications involve only the locations and brightnesses of the 
images, and thus require only the deflection and magnification components. Applications that 
involve the time delays, such as lensing measurements of the Hubble constant -f^o, also require the 
potential. 

If the mass distribution has circular symmetry, the deflection vector is purely radial and has 
an amplitude given by the 1-d integral 

1 M cyl (r) 



(j) r (r) 
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where M cy \(r) is the mass enclosed by a cylinder of radius r (the projected mass), which is often 
easily evaluated. The potential and magnification components can be obtained by integrating or 
differentiating 4> r . 

More general is the case of elliptical symmetry, in which the surface mass density has the form 

k = k(0, wh ere £ 2 = x 2 + y 2 /q 2 , (5) 



where q is the projected axis ratio and £ is an ellipse coordinate. This is the functional form in 
a coordinate system with the ellipse centered on the origin and aligned along the x-axis; other 
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coordinate systems can be reached by suitable translation and rotation. With elliptical symmetry 
the lensing properties can be written as a set of 1-d integrals (see Schramm 1990, although I have 
changed variables in the integrals), 

4>(x,y) = ^I(x,y) (6) 

4>x(x,y) = qxJ {x,y) (7) 

<l>y(x,y) = qyJi(x,y) (8) 

4>xx(x,y) = 2qx 2 K (x,y) + qJ (x,y) (9) 

4> yy {x,y) = 2qy 2 K 2 {x,y) +qJ 1 (x,y) (10) 

<Pxy(x,y) = 2qxyKi(x,y) (11) 



where the integrals are 



= /'iM *(«-)) (12) 
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J n (x,y) = [ — - ' — —TTrdu (13) 

Jo [l-(l-g>] n+1 / 2 V ' 

KJx,y) = / ' du 14 

^ ' Jo [1 - (1 - q 2 )u] n+1/2 

where £(u) 2 = u (x 2 + - _ J_ (15) 

and «'(£ 2 ) = dni^ 2 ) / d{^ 2 ) . Note from eq. (12) that the potential can be written as an integral over 
the circular deflection function <f> r from eq. (4), but 4> r must be evaluated at the appropriate ellipse 
coordinate £(u). 

All of the previous expressions assume that the surface density k is known. Some models 
have 3-d density distributions for which the projection integral cannot be evaluated analytically. 
In this case even a spherical lens model requires computationally expensive double integrals (the 
projection integral followed by the lensing integral). However, the double integrals can be rewritten 
as follows so the projection integral is replaced by the enclosed mass M(r), which can often be 
computed analytically. (The mass M{r) is the mass in spheres, which is different from M cy i(r) 
in eq. 4.) Writing re(r) as an integral over p(r) — the projection integral — and substituting 
p(r) = M'(r)/(47rr 2 ) where M'(r) = dM/dr, we find 

/s(r) = — — / du — v ; , (16) 
27rScr J r uVu 2 - r 2 

- ^h[ dv Thf K( r v^) +«'(^)] . <i7) 

where the second line represents a change of variables so the integral has a finite range, which is 
convenient for numerical integration. Combining eqs. (4) and (16) gives the circular deflection as 

1 f r f°° M'(v) 
<f>r = duu dv - ' = , (18) 



7i"S cr r J J u v\Jv 2 - u 2 
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where is the total mass. Integrating by parts then yields 

M(v) 
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(19) 
(20) 



(21) 
(22) 



where again the second line represents a change of variables for numerical integration. Eqs. (16) 
and (21), or eqs. (17) and (22) for numerical integration, represent the desired 1-d integrals for the 
surface density and deflection. The magnification also requires (p rr , which could be computed by 
differentiating (j) r ; however, it is easier to compute k and 4> r and then use the identity r _1 <ft r + (f) rr = 
2k to determine 4> rr . 



5. The Catalog 

Table 1 lists a wide variety of popular lens models, all of which are available in the gravlens 
software. This section summarizes what is known about the mass distributions and lensing prop- 
erties of the models. Analytic results are given where available, which includes all but one of the 
circular models and some of the elliptical models. Note that if the potential <p is regular at the 
origin, it is normalized to have <^>(0) = 0. Lensing is insensitive to an arbitrary constant added to 
the potential. 

I have tried to include relevant references. If you use results given here, please cite the original 
references rather than this catalog. As for the unreferenced results, some of them are new, while 
others are derived easily enough that references seem unnecessary. Use your own judgement about 
citing such material. 

Point mass: This model is inherently circular. A point mass M produces a lensing potential 

(f> = R 2 E \nr, (23) 
where the Einstein radius is (in angular units) 



Softened power law potential: Sometimes for analytic simplicity it is convenient to put 
the elliptical symmetry in the potential rather than in the density (e.g., Blandford & Kochanek 
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1987; Witt & Mao 1997, 2000; Evans & Hunter 2002), 

cj> = b(s 2 +x 2 + y 2 /q 2 ) a/2 -bs a . 

The deflection and magnification are given by simple derivatives. 

Softened power law ellipsoid: This model has a projected surface density 

1 b 2 - a 



(25) 
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(26) 



which represents a flat core with scale radius s, and then a power law decline with exponent a 
defined such that the mass grows as M cy \(r) oc r a asymptotically The core radius can be zero if 
a > 0. The model gives a softened isothermal model for a = 1, a modified Hubble model for a = 0, 
and a Plummer model for a = — 2 (see Binney & Tremaine 1987). The circular model has lensing 
properties 
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(27) 
(28) 
(29) 
(30) 



In the potential, 2F1 [o>, b; c; x] is a hypergeometric function, which can be written as or transformed 
into a quickly converging series (see Press et al. 1992; Gradshteyn & Ryzhik 1994, §9.1). Also, 
7£ = 0.577216... is Euler's constant, and ty(x) = d[lnT(x)]/dx is the digamma function, or the 
logarithmic derivative of the factorial function T(x). Analytic solutions for the elliptical model are 
possible for a = 0, ±1, and two of these are given below. Barkana (1998) gives a fast numerical 
algorithm for general softened power law ellipsoid models. 

Isothermal ellipsoid, a = 1: This model describes mass distributions with flat rotation 
curves (outside the core). Its lensing properties are: 



circular: 



(f>r = 

elliptical: (f> = 



r<j) r 
b 



bs In 
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(31) 
(32) 

(33) 
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; tan 



bq 



: tanh 



a/1 ~ q 2 x 

ip + s 



q 2 y 



ip + q 2 s 



(34) 
(35) 



where ip 2 = q 2 {s 2 + x 2 ) + y 2 . The elliptical solutions have been given by Kassiola & Kovner 
(1993), Kormann, Schneider &: Bartelmann (1994), and in the simple form quoted here by Keeton 
& Kochanek (1998). In the limit of a singular (s = 0) and spherical (q = 1) model, b is the Einstein 
radius of the model and is related to the 1-d velocity dispersion a by 
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(36) 



(in angular units). 

a = — 1 ellipsoid: This model corresponds to an unnamed density profile with £ oc r~ 3 
(p oc r~ 4 ) asymptotically. Its lensing properties are: 



circular: 



In 



_|_ g 2 _|_ r 2 

2s 



0r = 



1 - 



\/s 2 + r 2 



elliptical: 
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where ip 2 = q 2 (s 2 + x 2 ) + y 2 . The elliptical solutions are given by Keeton & Kochanek (1998). 



(37) 

(38) 

(39) 
(40) 
(41) 



Pseudo-Jaffe ellipsoid: A standard Jaffe (1983) model has a 3-d density distribution p oc 
r~ 2 (r + a)~ 2 where a is the break radius. For lensing it is useful to modify this model and write 
p oc (r 2 + s 2 )~ l {r 2 + a 2 ) -1 , where a is again the break radius and we have added a core radius 
s < a. The projected surface density of the elliptical model has the form 



(42) 



which is constant inside s, falls as R^ 1 between s and a, and falls as i?~ 3 outside a; the total mass 
is M = irY, cr qb(a — s). Eq. (42) defines the pseudo- Jaffe ellipsoid. In the limit a — > oo it reduces 
to the isothermal ellipsoid (a = 1). In the limit a — ► s it reduces to the a = —I ellipsoid, although 
the limit must be taken in p rather than in k (i.e., the limit must be taken before the projection 
integral is evaluated). The pseudo- Jaffe model is equivalent to a combination of two softened 
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isothermal ellipsoids, so its lensing properties can be computed with appropriate combinations of 
eqs. (33)-(35). 

King model: The King model can be approximated as a combination of two softened isother- 
mal models (see Young et al. 1980; Barkana et al. 1999), 

2.126 1.756 



«(0 



V0.75r2 + £2 ^2.997-2 + ^ 



(43) 



It has a single scale radius r s . This approximation is convenient because it is written as the 
difference of two softened isothermal ellipsoids, so its lensing properties can be computed with 
appropriate combinations of eqs. (33)-(35). 

De Vaucouleurs model: This is the prototypical constant mass-to-light ratio lens model (de 
Vaucouleurs 1948), with surface mass density 



k(0 = K exp \-k^/R e )^ 4 



(44) 



where k = 7.66925001 and R e is the major-axis effective (or half-mass) radius. The circular 
deflection is (Maoz & Rix 1993) 
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where ( = k (r/i? e ) 1//4 . The elliptical model can be computed numerically with eqs. (6)— (11). 

Hernquist model: The Hernquist (1990) model is a 3-d density distribution with a projected 
distribution that mimics the luminosity distribution of early-type galaxies. It has the form 



Ps 



(r/r s )(l + r/r s ) 3 



(46) 



where r s is a scale length and p s is a characteristic density. The projected surface mass density has 
the form 



-3 + (2 + x 2 )F(x)} 



(x 2 - l) 2 

where x = r/r s , k s = p s r s /S cr , and J r (x) is the function 

' -, I , tan~Va; 2 - 1 (x > 1) 
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(48) 



A useful technical result is the derivative of this function, 

1 - x 2 T{x) 
x(x 2 — 1) 
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(49) 
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The circular potential and deflection are 



<f> = n s r s 



x 



4>r 



2 K/a V , 



x[l - F{x)} 
x 2 -l ' 



(50) 
(51) 



where again x = r/r s . The elliptical model can be computed numerically with eqs. (6)— (11). 

NFW model: Cosmological iV-body simulations (e.g., Navarro et al. 1996, 1997) suggest that 
dark matter halos can be described by a "universal" density profile with the form 

(52) 



(r/r s )(l + r/r s ) 2 



For the spherical NFW model, the projected surface mass density, potential, and deflection arc 
(Bartelmann 1996; Meneghetti, Bartelmann k Moscardini 2001, 2002; Golse & Kneib 2001) 
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(53) 
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(55) 



where x = r/r s , k s = p s r s /S cr , and the function F{x) is the same as in the Hernquist model. The 
elliptical model can be computed numerically with eqs. (6)— (11). 

Cuspy NFW model: Moore et al. (1998, 1999) have suggested that the inner cusp of the 
NFW profile is too shallow, so Jing & Suto (2000), Keeton & Madau (2001), and Wyithe et al. 
(2001) have studied a generalized NFW-type profile of the form 



Ps 



(r/r s p(l+r/r s f-~f ' 



(56) 



so the central cusp has p oc r~ 7 . The projected surface density cannot be computed analytically 
even for a spherical halo. For a spherical model, eqs. (16) and (21) allow the surface density and 
deflection to be written as 
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I -^ 2 F 1 [3-7,3- 7 ;4-7;-^]+ f'dyiy + xy- 3 - ^ ^ 1, 
[3-7 J y J 

where x = r/r s , k s = p s r s /T, cr , and 2^1 is the hypergeometric function. The integrals can be 
evaluated for several values of 7 to obtain simple analytic expressions. For 7 = 1 the model reduces 
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to the standard NFW model (see eqs. 53 and 55). Two other analytic cases are as follows: 

7 = 2 : «(r) = k s ^ - 2JF(x) , (59) 

0r = 4 Ks r s \ + \ In (|) +i^^(x) , (60) 
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Elliptical NFW-cusp models can be computed numerically with eqs. (6)— (11). In general 
they are slow to compute because they require double integrals. 

Cuspy halo models: To obtain a general cuspy model that is more amenable to lensing, 
Muhoz et al. (2001) introduce a model with a profile of the form 

P = ( r /r s )7[l + ( r /r s ) 2 ]( n -T)/ 2 ' ^ 

where again r s is a scale length, and 7 and n are the logarithmic slopes at small and large radii, 
respectively. This model is a subset of the models whose physical properties were studied by Zhao 
(1996). The central cusp must have 7 < 3 for the mass to be finite. For (7, n) = (1,4) this is a 
pseudo-Hernquist model, for (1,3) it is a pseudo-NFW model, and for (2,4) it is a singular pseudo- 
Jaffe model. Compared with eq. (56), replacing (1 + r/r s ) with + (r/r s ) 2 does not greatly 
change the profile shape but does make it possible to solve the spherical model analytically (Muhoz 
et al. 2001), 
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(64) 
(65) 



2 ' 2 ' 2 ' 1 + x 2 



where x = r/r s , k s = /? s r s /S cr , 2-^1 is the hypergeometric function, B(a, b) = T(a)T(b)/T(a + 6) is 
the Euler beta function. For n = 3 eq. (65) is not valid (it has singularities of the form 0/0), so it 
must be replaced with one of the two following expressions: 
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where is the digamma function (see eq. 27), and G(b, c; z) is the function 

^ 2 F 1 [a,b,c; z] - 1 
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(68) 
(69) 

(70) 



Note that eqs. (66) and (67) are equivalent and can be transformed into each other using identities 
for the hypergeometric function (see Gradshteyn & Rhyzik 1994, §9.131). It is convenient to use 
eq. (66) for x < 1 and eq. (67) for x > 1 in order to make the series for the G function converge 
rapidly. 

Elliptical cusp models n(£) can be computed numerically with eqs. (6)-(ll), or with a Fourier 
series solution (Chae 2002). 

Nuker law: Many early-type galaxies have surface brightness profiles that can be modeled as 
a Nuker law (e.g., Lauer et al. 1995; Byun et al. 1996), 





-7 




(k) 













( 7 -/3)/a 



(71) 



where 7 and f3 are inner and outer power law indices, respectively, is the radius where the break 
in the power law occurs, a gives the sharpness of the break, and lb is the surface brightness at the 
break radius. If the luminosity distribution has circular symmetry and the mass-to-light ratio is T, 
the lensing deflection is (Keeton 2002) 



2 l+(/3-7)/a 

2-7 



KbH — 

Jb 



1-7 



2-^1 



a 



7 0-7. 1 + 2 



7. 



a 



a 



(72) 



where Kb = TIb/T, CT is the surface mass density at the break radius in units of the critical density for 
lensing, and 2 -Pi is the hypergeometric function. Ellipsoidal models can be computed numerically 
with eqs. (6)-(ll). 



Exponential disk: The projected surface density is 

«(0 = (f 1 K o exp [-£/Rd] , 



(73) 



which represents a thin, circular disk with intrinsic central density kq and scale length R^, seen in 
projection with axis ratio q = \ cosi\ where i is the inclination angle (such that i = 0° is face-on 
and i = 90° is edge-on). The circular deflection is 



p 2 r 

r 



1 



1 



— ] e- r / Rd 

Rd 



(74) 
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The elliptical model k(£) can be computed with eqs. (6)-(ll), or it can be approximated with one 
or more Kuzmin disks (see Keeton & Kochanek 1998). 

Kuzmin disk: The a = — 1 ellipsoid can be re-interpreted as the projection of a thin disk, 
in which case it corresponds to a Kuzmin (1956) or Toomre (1962) Model I disk; see Keeton Sz 
Kochanek (1998). Its projected surface density is 

^)=^ 1 -or s 3 (r, 2 + e 2 )- 3/2 , (75) 

where is the intrinsic central surface density of the disk, and q = | cos i\ is again the projected 
axis ratio of the inclined disk. The only difference between the Kuzmin disk and the a = — 1 
ellipsoid is the normalization. 

External perturbations: Objects near the main lens galaxy or along the line of sight often 
perturb the lensing potential. If the perturbation is weak it may be sufficient to expand the 
perturbing potential as a Taylor series and keep only a few terms. In a coordinate system centered 
on the lens galaxy, the expansion to 3rd order can be written as (see Kochanek 1991; Bernstein Sz 
Fischer 1999) 



r 



r 



3 



<p + b • x + — k - 7 cos 2(0 - 7 ) + — 5 cos(6 - S ) - e cos 3(0 - e ) + . . . (76) 



The Oth order term </>o represents an unobservable zero point of the potential and can be dropped. 
The 1st order term b • x represents an unobservable uniform deflection and can also be dropped. 
The 2nd order term k represents the convergence from the perturbing mass and is equivalent to a 
uniform mass sheet with density S/S cr = k. The only observable effect of this term is to rescale 
the time delay(s) by 1 — k, which leads to the "mass sheet degeneracy" (e.g., Falco, Gorenstein 
& Shapiro 1985); hence this term is often omitted from lens models and introduced a posteriori 
using independent mass constraints (see, e.g., Bernstein & Fischer 1999). The 2nd order term 7 
represents an external tidal shear with strength 7 and direction 7 . The 3rd order term 5 arises 
from the gradient of the surface density k(x) of the perturber; it has an amplitude 5 = (3/4)|Vre| 
and a direction equal to the direction of V/c. The 3rd order term e arises from the m = 3 multipole 
moment of the perturbing mass. The constant coefficients (k, 7, S, e) are all evaluated at the position 
of the lens galaxy, and the corresponding direction angles are written here as theory angles measured 
counter-clockwise from the x-axis. 
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Table 1. Mass Models for Lensing 





N 


lion oi "i"^ r /~» f T* 1 




Point mass 


o 


<5(x) 


<5(x) 


Power law or a-models 


2 


( s 2 +r2) (-3)/2 


( s 2 + r2) (-2)/2 


Isothermal (a = 1) 


1 


(s 2 +r 2 ) _1 


( s 2 + r 2)-V2 


a = — 1 


1 


(s 2 + r 2 )- 2 


( s 2 + r 2)"3/2 


Pseudo-Jaffe 


2 


( S 2 +r 2)-l( a 2 +r 2)-l 


( s 2 + r 2)-V2_ ( a 2 +r 2)"l/2 


King (approximate) 


1 




2.12(0.75r 2 + r 2 )" 1/2 
-1.75 (2.99r s 2 + r 2 )" 1/2 


de Vaucouleurs 


1 




exp [-7.67(r/ J R e ) 1 / 4 ] 


Hernquist 


1 


r -1 (r s + r)~ 3 


see eq. (47) 


NFW 


1 


r _1 (r s + r)~ 2 


see eq. (53) 


Cuspy NFW 


2 


r~ 7 (r s + r) 7_3 


see eq. (57) 


Cusp 


3 


r -7 ( r 2 +r 2)(7-n)/2 


see eq. (64) 


Nuker 


4 




see eq. (71) 


Exponential disk 


1 




exp[-r/R d ] 


Kuzmin disk 


1 




(r s 2 + r 2 )- 3/2 



Note. Density profiles for lensing mass models; see §5 for detailed definitions, 

including normalizations. Three-dimensional density profiles are not given for the King, 
de Vaucouleurs, Nuker, exponential disk, and Kuzmin disk models because these models 
are defined by their surface densities. Column 2 indicates the number of parameters 
associated with the radial profile alone; each model would also have parameters for the 
position and the mass scale, and elliptical models would have parameters for the ellipticity 
and orientation. The profiles are given for spherical models; elliptical models are defined 
by k(£) where £ is an ellipse coordinate (see eq. 5). 



